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Comparative tests of isospin-symmetry-breaking corrections to superallowed + — >0 + 

nuclear j3 decay 

I.S. TowneiQ and J.C. Hard}fl 

Cyclotron Institute, Texas A&M University, College Station, Texas 77843 
(Dated: November 29, 2010) 

We present a test with which to evaluate the calculated isospin-symmetry-breaking corrections 
to superallowed + —*■ + nuclear /3 decay. The test is based on the corrected experimental Tt 
values being required to satisfy conservation of the vector current (CVC). When applied to six sets 
of published calculations, the test demonstrates quantitatively that only one set - the one based 
on the shell model with Saxon- Woods radial wave functions - provides satisfactory agreement with 
CVC. This test can easily be applied to any sets of calculated correction terms that are produced 
in future. 

PACS numbers: 23.40.Bw, 23.40.Hc 



I. INTRODUCTION 

Superallowed + — >• + (3 decay between T = 1 nuclear 
analog states has been a subject of continuous and of- 
ten intense study for six decades. The ft values for such 
transitions are nearly independent of nuclear-structure 
ambiguities and depend uniquely on the vector part of 
the weak interaction. Their measurement gives us access 
to clean tests of some of the fundamental precepts of 
weak-interaction theory, and, over the years, this strong 
motivation has led to very high precision being achieved 
both in the experiments and in the theory used to inter- 
pret them. 

The most recent survey of world data [l| finds ten of 
these superallowed transitions with measured ft values 
known to 0.1% precision or better, and three more that 
have a precision of between 0.1 and 0.3%. An analysis 
of the ft values [lj demonstrated that the vector cou- 
pling constant, G v , has the same value for all thirteen 
transitions to within ±0.013%, thus confirming a key 
part of the Conserved Vector Current (CVC) hypothesis; 
and it set an upper limit on a possible scalar current at 
0.2% of the vector current. With both these outcomes 
established, the results could then be used to extract 
a value for V u( i, the up-down element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, with which the top- 
row unitarity test of the CKM matrix yielded the result 
Q 0.9999(6). This is in remarkable agreement with the 
Standard Model, and the tight uncertainty significantly 
limits the scope for any new physics beyond the model. 
Further tightening of the uncertainty would increase the 
impact of this result even more. 

Although the role played by nuclear structure is rel- 
atively small, the precision currently reached by experi- 
ment is such that the theoretical uncertainties introduced 
by correction terms required in the analysis of the ft- 
value data now predominate over the experimental un- 
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certainties. Two of these correction terms depend on nu- 
clear structure and together they are the second largest 
contributor to the overall uncertainty in V uc i- The largest 
contributor is the nucleus-independent component of the 
radiative correction but at present there seems little op- 
portunity for further improvement there. 

Thus it is the nuclear-structure dependent terms that 
have attracted the greatest attention, particularly re- 
cently. The most widely used of these latter correction 
terms are those calculated by the present authors, which 
have been tabulated for all the superallowed transitions 
of interest in Ref . [3J . However, there are a growing num- 
ber of alternative choices [J-Q available for one of the 
two correction terms - the one that accounts for isospin 
symmetry breaking - including a set we offer ourselves 
[l|. There has also been a claim, albeit unsupported by 
any detailed computations, that our calculations neglect 
a radial excitation term, which is purported to be im- 
portant (lfj. To counterbalance that, however, there are 
two recent papers that confirm our result: one [ll| does 
so based on a semi-empirical analysis of the data, while 
the other [12| quotes the average results from a Skyrme- 
density-functional-theory calculation in which simultane- 
ous isospin and angular-momentum projection has been 
incorporated. 

Clearly it would be valuable if the various sets of calcu- 
lated isospin-symmetry-breaking correction terms could 
be tested against the data, and their relative merits quan- 
titatively evaluated, since this must surely be a first step 
in any attempt to reduce the uncertainty attributed to 
these corrections. In this paper, we address ourselves to 
devising and then applying such a test. 

We begin by describing how information on the fun- 
damental weak-interaction parameters is extracted from 
the experimental ft- value data. We will overview the role 
played by all the theoretical corrections but will focus, in 
particular, on the isospin-symmetry-breaking term. This 
will lead naturally to the test we propose as a means 
of evaluating the efficacy of any calculated set of these 
terms available now or in the future. We will then out- 
line the methods currently used to calculate the isospin- 
symmetry-breaking term, and proceed to apply our test 



to each of them. Finally we will evaluate the results of 
the test and present our conclusions. 



II. THE ANALYSIS OF SUPERALLOWED 
BETA TRANSITIONS 

Supcrallowed Fermi beta decay between + states de- 
pends uniquely on the vector part of the hadronic weak 
interaction. According to CVC, when the decay occurs 
between isospin T = 1 analog states the measured ft val- 
ues should be the same irrespective of the nucleus, viz. 



ft 



K 



Gl\M F \ 



const, 



(1) 



where K/{hcf = 2TT 3 h\n2/(m e c 2 ) 5 = (8120.2787 ± 
0.0011) x 10 -10 GcV~ 4 s; G v is the vector coupling con- 
stant for semi-leptonic weak interactions; and Mp is the 
Fermi matrix element. The CVC hypothesis asserts that 
the vector coupling constant, G v , is a true constant and 
not renormalised to another value in the nuclear medium. 
In practice, Eq. (TTJ) has to be amended slightly. Firstly, 
there are radiative corrections because, for example, the 
emitted electron may emit a bremsstrahlung photon that 
goes undetected in the experiment. Secondly, isospin is 
not an exact symmetry in nuclei so the nuclear matrix 
element, Mp, is not the same for all superallowed tran- 
sitions but is slightly reduced from its ideal value by a 
different amount in each case. This leads us to write: 



\Mp 



\M \ 2 (1-6 C ), 



(2) 



where Mq is the exact-symmetry value, which for T = 
1 states is Mq = \/2; and 8c is the isospin-symmetry- 
breaking correction, which takes on a different (small) 
value for each transition. Thus, we define a "corrected" 
Tt value as 



Tt = ft(l + 8 R )(l + 8 NS -8c) = 



K 



2G2(1 + A X ) : 



(3) 



where, in addition to the terms already defined, A^ is 
the transition-independent part of the radiative correc- 
tion; while the terms 8' R and 8ns comprise the transition- 
dependent part of the radiative correction, the former be- 
ing a function only of the electron's energy and the Z of 
the daughter nucleus, while the latter, like 8c, depends 
in its evaluation on the details of nuclear structure. 

From this equation, it can be seen that a single mea- 
sured transition establishes a value for Tt, and hence 
G v . This result could, in principle, then be used to de- 
termine Vud via the relationship V u d = G V /G F , where G F 
is the well known weak-interaction constant for muon de- 
cay [13| . However, a value for V u d derived from a single 
superallowed transition would be reliant upon a single 
pair of structure-dependent correction terms, 8ns an d 
8c, without there being any independent verification of 
those terms' validity; so, in practice, as many transitions 
as possible arc measured and their resultant Tt values 



compared. If they satisfy CVC by being statistically con- 
sistent with each another, then one is justified in taking 
an average value of J-t, from which G v and V u d can then 
be derived. 

If they arc not consistent with each other, then one 
can proceed no further since inconsistency must signal a 
failure either of the calculated structure-dependent cor- 
rections or else of the CVC hypothesis itself. In cither 
case, an average value of Tt has no defined significance 
and certainly cannot be used to obtain a value for V u d- 

Here we find the basis for a test of the calculated 
structure-dependent correction terms: How well do they 
do in producing a consistent set of Tt values from the ex- 
perimental ft values? The latter show very pronounced 
differences from one transition to another, and the extent 
to which those differences are successfully removed by a 
given set of calculated correction terms would be a sensi- 
tive measure of the efficacy of the calculations involved. 
Naturally, such a test is only as good as the CVC hypoth- 
esis. However, we believe that most would agree that a 
persistent scatter in the derived Tt values is more likely 
to be due to a deficiency in the calculated corrections 
rather than to a failure of CVC. 



III. THE TEST 

Our test is based upon the premise that CVC is valid at 
least to ±0.03%, which is the level of precision currently 
attained by the best /t-value measurements. Under that 
condition, a valid set of structure-dependent correction 
terms should produce a statistically consistent set of Tt 
values, the average of which we can write as Tt. It then 
follows from Eq. ([3]) that, for each individual transition 
in the set, we can write 



Tt 



ms 



ft(i + s' R y 



(4) 



For any set of corrections to be acceptable, the calcu- 
lated value of 8c ~ 8ns for each supcrallowed transition 
must satisfy this equation, where ft is the measured re- 
sult for that transition and Tt has the same value for all 
of them. Thus, to test a set of correction terms for n 
superallowed transitions, one can treat Tt as a single ad- 
justable parameter and use it to bring the n results from 
the right side of Eq. @ , which are based predominantly 
on experiment, into the best possible agreement with the 
corresponding n calculated values for 8c — 8ns- The nor- 
malized x 2 , minimized by this process, then provides a 
figure of merit for that set of calculations. 

As it happens, there is only one set of calculations 
available for Sns @: 13 but many for the isospin- 
symmetry-brcaking term 8c- It therefore becomes more 
useful to rearrange Eq. (f4]) to read: 



Tt 



Jc 



?NS 



ft{i + s' R y 



(5) 



The same least-squares minimization process can of 
course be used in the application of this equation. 



IV. AVAILABLE CALCULATIONS FOR S c 

There have been a number of methods used over the 
years to calculate the isospin-symmetry-breaking correc- 
tion to superallowed /3 decay. We describe some of them 
here, in chronological order. 



A. Damgaard model 



which, for the light nuclei we are interested in, exhibits 
the general behaviour 5 C oc A 4 / 3 with some shell struc- 
ture superimposed through the choice of oscillator quan- 
tum numbers n and L In particular, a proton radial 
function with one radial node gets a factor of two en- 
hancement in its 8c value over one that has no radial 
nodes simply from the factor (n + 1) in Eq. (|11[) . 

We have used Eq. (fTTj) to derive 5c values for the thir- 
teen best known superallowed transitions. These transi- 
tions are listed by parent nucleus in the first column of 
Table HI and the 5c results for this model appear in the 
fifth column of the same table. We will use these results 
in our comparative tests of all models. 



The first model was proposed in 1969 by Damgaard 
[J] and was improved 8 years later by Towner, Hardy 
and Harvey [nf . The idea is that the proton involved in 
beta decay has a different radial wave function than the 
neutron into which it transforms because it is influenced 
by the Coulomb interaction with all the other protons in 
the nucleus. If the other protons present a uniform charge 
distribution of radius, R, then the Coulomb interaction 
for a proton at r < R is 

7p 2 A 
V c (r) = ~ J2 (H - W) (I - **») S(r-n), (6) 
i=l 

where t z {i) = — i if nucleoli i is a proton, and = +| if it 
is a neutron. 

Using an oscillator model as a basis, Damgaard ex- 
panded the proton radial function in terms of a complete 
set of neutron oscillator functions. The set comprised 
states of the same orbital angular momentum, t, but dif- 
fering numbers of radial nodes, n. Most of the mixing 
turned out to be with the state with one more radial 
node, so 



u pro on^ ^ ^ _ ^2)1/2^ ^ r ) + mi„ +M (r). 



(7) 



The mixing amplitude comes from first-order perturba- 
tion theory, 



a = {u n +i,t\V c \un t i)(AE) 1 

AE = E n+ i_( — E„ p = 2tku, 



(8) 



and the Fermi matrix element between T = 1 states is 
given by 



\M F \ 2 = 2(l-a 2 ). 



(9) 



Upon evaluating the Coulomb matrix element using os- 
cillator functions with V c taken from Eq. ([5]), Damgaard 
obtained 



5 C = a 2 = 



4*4 



e 4 h 



{huj) 4 R 6 16m 2 



(n + l)(n + £+§). (10) 



If we adopt the relationships Tuo = 41A 1 / 3 MeV and 
R = 1.2 A 1 / 3 fm, this expression becomes 



5 C = 0.2645 x Z 2 A~ 2/3 (n + l)(n + 1 + §), 



(11) 



B. Shell model with Saxon- Woods radial wave 
functions (SM-SW) 

This model was introduced by Towner, Hardy and Har- 
vey in 1977 [15| and improved upon several times since 
then |3j, [l6[ . In their approach, the Fermi matrix element 
is defined by 



Mr 



E 



(fW a K\ 



E 



(/|at»W& Q |i), (12) 



where a^ creates a neutron and b a annhilates a proton in 
state a. Here \i) and |/) are the exact ^4-body state vec- 
tors for the full Hamiltonian, and \tt) represents a com- 
plete set of (A — l)-body parent states. If this Hamilto- 
nian commutes with the isospin operators, then \i) and 
|/) are exact isospin analogues of each other, and the 
symmetry-limit matrix element is 



M = J2\(M 



(13) 



which for T = 1 states corresponds to Mo = v2. How- 
ever, with isospin not being an exact symmetry, \i) and 
|/) are not exact isospin analogues; nevertheless the re- 
sulting matrix element Mp is not very different from Mo , 
the relationship between them being given by Eq. (|2|): 
viz. M F = Mq(1 — 5c), where 5c is small. 

Ideally, to obtain 5c one would compute Eq. f|L2[) us- 
ing the shell model, and introduce Coulomb and other 
charge-dependent terms into the shell-model Hamilto- 
nian. However, the shell-model space would have to be 
huge to include all the potential states with which the 
Coulomb interaction might potentially connect. Since 
this is not a practical proposition, a model approach was 
developed in which 5c is divided into two parts: 



$c = Sci + <5(72- 



(14) 



For Jci, one computes 

Y,(7W a \n)(n\b a \J) = M (l - <fci) 1/2 , (15) 



where |z) and |/) are not the exact eigenstates that 
appear in Eq. (fT2|) but are the shell-model eigenstates 
of an effective Hamiltonian (including charge-dependent 
terms) evaluated in a tractable shell-model space. How- 
ever, this space is not large enough to allow for mixing 
with functions having a different number of radial nodes, 
so the term 8c2 is introduced to compensate for that 
limitation. This second term is derived from 



£l(/kk)l 2 ^ = M (l - Sc%) 



1/2 



(16) 



where r^ is a radial overlap integral of proton and neu- 
tron radial functions. If the proton and neutron radial 
functions were identical, then it would follow that r* = 1, 
and (5(72 = 0. But, since they are not identical, a finite 
correction 8ci is obtained. The idea is that nodal mixing 
mainly impacts on the radial functions - as demonstrated 
by the Damgaard model - and so its impact is best mod- 
elled by Eq. (Unj)- 

The wave function for the decaying A-body state, \i), 
is expanded in a set of parent states of (A — l)-nucleons, 
\tt), plus a proton; while that of the daughter A-body 
state, |/), is expanded in terms of the same set of parent 
states plus a neutron. The expansion coefficients are ob- 
tained from a shell-model calculation. Isospin-symmctry 
breaking is introduced by allowing the radial function for 
the proton in these expansions to differ from that of the 
neutron. In this model, these radial functions are taken 
to be cigenfunctions of a Saxon- Woods potential. The 
well-depths of the proton and neutron potentials are ad- 
justed so that the asymptotic forms of the radial function 
go as e~ ar , where a 2 = 2mS/h 2 , with m being the nu- 
cleon mass, and S the experimental separation energy 
for the proton (or neutron) in the A-body state. Further 
details can be found in Q. 

It is important to realize that this model is really semi- 
phenomcnological in its application. In addition to the 
match with experimental separation energies in the cal- 
culation of (5c2, the radius of the Saxon- Woods potential 
in each case was set to the value determined experimen- 
tally for the charge radius by electron scattering |16| , and 
the shell-model parentage was linked to measured single- 
nucleon transfer reactions Q . The value of Sci was also 
constrained by comparison with experiment. First, for 
each superallowed transition the single-particle energies 
of the proton orbits were shifted relative to the neutrons, 
the exact amount being determined from the spectrum of 
single-particle states in the closcd-shell-plus-proton ver- 
sus the closcd-shell-plus-neutron nucleus. Second, the 
two-body Coulomb interaction among the valence pro- 
tons was adjusted in strength for each decay so that 
the measured b coefficient in the isobaric multiplet mass 
equation (IMME) was exactly reproduced for the multi- 
plet involved in that decay. Third, the charge-dependent 
nuclear interaction, which had been incorporated by a 
^2% increase in all the T = 1 proton-neutron matrix el- 
ements relative to the neutron-neutron ones, was tuned 
to give agreement with the measured c coefficient of the 



IMME. 

The current best values for 8c as calculated with this 
model are listed in the fifth column of Table VI in Ref. [3| , 
and arc reproduced here in the sixth column of Table [I] 



C. Shell model with Hartree-Fock radial wave 
functions (SM-HF) 

Beginning in 1985, Ormand and Brown [5|, Iff adopted 
the same general procedure as the one just described, 
splitting (5c into two components, the first of which, 8ci, 
incorporated configuration mixing within a restricted 
shell- model space, and the second, 8c2, accounted for 
mixing with all other states by evaluating the mismatch 
in the parent and daughter radial wave functions. The 
shell-model aspects of their model were the same as the 
SM-SW model, but their radial functions were taken to 
be cigenfunctions of a mean-field Hartree-Fock potential 
rather than of a Saxon- Woods potential. As in the SM- 
SW model, the strength of this mean field was readjusted 
so that the asymptotic forms of the proton and neutron 
radial functions were matched to their respective separa- 
tion energies. 

Ormand and Brown's [5j, [y] protocol was to perform 
two Hartree-Fock calculations with a Skyrme interaction: 
one for the decaying A-body state, whose mean field pro- 
vided the proton function, and the other for the daugh- 
ter A-body state, whose mean field provided the neutron 
function. However, it was noted more recently by Hardy 
and Towner [1| that there is a problem with this protocol: 
the Coulomb part of the proton mean field has asymp- 
totically the wrong form, falling off as (Z + l)e 2 /r rather 
than Ze 2 /r. They therefore modified the protocol to 
just a single Skyrme-Hartree-Fock calculation performed 
in the (A — l)-body state, whose mean field provided for 
both the proton and neutron radial functions. In this pro- 
cedure, the Coulomb interaction automatically has the 
right asymptotic form. Further details can be found in 

We have obtained the 8c values for this model by 
adding the "adopted" 8ci numbers from Table III of 
Ref. [3[ (the same as we used for the SM-SW model) 
and the "HF" (5 C2 numbers from Table XI of Ref. 0). 
The results appear in column seven of our Table [T] 



D. Hartree-Fock with Random Phase 
Approximation (RHF-RPA and RH-RPA) 

In 1996, Sagawa, Van Giai and Suzuki [7| introduced 
a new model, in which a Skyrme-Hartree-Fock calcu- 
lation was performed for each even-even A-body sys- 
tem: the parent for the cases in which the superallowed 
decay proceeds from a T z = —1 parent nucleus, and 
the daughter for cases of decay from a T z = parent 
nucleus. The odd-odd nucleus was then treated as a 
particle-hole excitation built on the even-even Hartrce- 



Fock state. The particle-hole calculation was carried 
out in the charge-exchange random-phase approximation 
(RPA) in a model space extending up to lOfuv excita- 
tion, with radial functions up to five nodes. The low- 
est state in the RPA spectrum was identified as the 
isobaric analogue state - the state actually involved in 
the superallowed Fermi beta decay. Unlike the previous 
two methods, there was no adjustment to reproduce ex- 
actly the energy of the analog state, but the authors did 
check that their results were typically within 500 keV 
of the experimental value. Isospin-symmetry breaking 
was introduced by the presence of a Coulomb interac- 
tion, augmented by explicit charge-symmetry-breaking 
and charge-independence-breaking interactions included 
in the two-body force used in the Hartree-Fock calcula- 
tion. 

Since first results from this model appeared [7| signifi- 
cant progress has been made in self-consistent RPA calcu- 
lations in charge-exchange channels. Skyrme zero-range 
interactions have been replaced by finite-range meson- 
exchange potentials involving a, u), p, and ir mesons, and 
a relativistic rather than a non-relativistic treatment can 
be used. In 2009 Liang, van Giai and Mcng @ pub- 
lished improved results from relativistic Hartree-Fock, 
RHF-RPA, calculations with three different effective in- 
teractions, as well as from relativistic Hartrcc (only), 
RH-RPA, calculations with density-dependent meson- 
nucleon couplings and non-local interactions. The results 
were not particularly sensitive to the interaction used, so 
in performing our tests, we just use one interaction for 
each type of calculation: PKOl for RHF-RPA and DD- 
ME2 for RH-RPA §. The corresponding values for S c , 
were taken from Table I in Ref. [8|, and are reproduced 
here in columns eight and nine of Table |U Note that the 
authors of Ref. [8[ only calculated Sc values for 8 of the 
13 well known superallowed transitions. 



E. Isovector Monopole Resonance (IVMR) 

In 2009, Auerbach Q introduced a model in which he 
assumed that isospin-symmetry breaking in superallowed 
/3 decay is due entirely to mixing with the giant monopole 
state. 

The isovector part of the Coulomb interaction, which 
appears in Eq. ((5]) , is defined to be the isovector monopole 



operator Mg , thus 



M«=5>^), 



(17) 



r(l) 



where A/q is a spherical tensor in isospin space of rank 
1, with its z-component equal to 0. We write \M) to 
be the giant monopole state, which is created by the 
application of operator Mg to the ground state. If 
the ground state has N=Z with isospin quantum num- 
bers T=0 and T z =0, then the giant monopole state is 
a unique state with quantum numbers T= 1 and T z = 0. 



But if the ground state has a neutron excess, with 
T=T Z = |(iV — Z), then the monopole state is split into 
two components, one with isospin T and the other with 
isospin T+l. In this case the ground-state wave function 
is designated by \T,T) and the two components of the 
monopole state by \Mt,t) an d \Mt+i,t)- Furthermore, 
the isobaric analogue of this ground state, \T, T— 1), has 
its giant monopole state split into three isospin compo- 
nents, \Mt-i,t-i), \Mt,t-i) and \Mt+i,t-i)- 

By assuming that the giant monopole state is the sole 
source of isospin-symmetry breaking in superallowed de- 
cays, Auerbach [9| could write the wave functions for the 
two states involved in the /3 decay as 



l*i) 

l*2> 



(|r,T) + eo|MT,T) + e 1 |Mr+i 1 r))wr 1 

(|T,T-l)+»/_ 1 |AfT_i 1 r-i) 
+ffo|M r ,r-i) +7 ? i|Mr+i,T-i})iV 2 - 1 , (18) 



where 



N 2 
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rii) 



(19) 



and amplitudes e^ and rn can be expressed via pertur- 
bation theory in terms of Coulomb matrix elements be- 
tween the ground state and the respective components 
of the isovector monopole state. Based on this result, 
he then derived the corresponding isospin-symmetry- 
breaking correction to superallowed /3 decay, which he 
wrote to order 0(e 2 ,f] 2 ) as 



fc = V-i + ^o-Vof 



-e 1 +?7 1 -2ei?7i 



2T+1 
T 



1/2 



(20) 



Auerbach next argued that the Coulomb matrix ele- 
ments of differing isospins are all related to each other 
via isospin Clebsch-Gordan coefficients. He thus found 
that the coefficients eo, 7?_i, r]o and rji could all be ex- 
pressed in terms of the one isospin-mixing amplitude t\. 
In this way, the expression for Sc, Eq. (|20)) . reduces to 



S C = 8(T+1) 



r, 



£%ujA 



(21) 



where £ is related to the particle-hole interaction energy 
required to place the centroid of the giant monopole res- 
onance at the appropriate energy; and V\ is related to 
the strength of the symmetry potential that sets the en- 
ergy splitting between the components of the monopole 
state. Auerbach chose V\ = 100 MeV, £ = 3, and 
fibj = 41A -1 / 3 MeV and estimated e\ by appealing to a 
number of "gross" models discussed in Ref. [17| : a hydro- 
dynamical model, models based on non-energy-weightcd 
and energy-weighted sum rules, and a microscopic model. 
Each enabled him to obtain a simple expression for Sc 



TABLE I: Input data for the tests of the isospin-symmetry-breaking corrections, 5c, obtained from the various models described 
in Sect. lIVl The experimental ft values come from Table IX in the most recent survey of world data [J; however, in order to 
ensure undiluted normal statistics, we have set all the "scale factors" used in that reference equal to 1, with the consequence 
that the uncertainties quoted for most cases are smaller than those listed in Ref. [l[. The calculated values of 5' R and 5ns 
come from Table VII of Ref. [3|. The 5c values tabulated in the last six columns were obtained as follows: the "Damgaard" 
values were derived from our Eq. (|11[) ; those labeled SM-SW came from Table VII of Ref. [3|; the SM-HF values were obtained 
by adding the "adopted" 5ci numbers from Table III of Ref. [3] and the "HF" 5c2 numbers from Table XI of Ref. [J|; the 
RHF-RPA values were taken for the PKOl effective interaction given in Table I of Ref. [g|; the RH-RPA numbers, which 
corresponded to the density-dependent DD-ME2 effective interaction, were taken from the same table and reference; the IVMR 
values were calculated from our Eq. (|22[) . which is the same as Eq. (32) in Ref. [S|. 



Parent 


Experimental 
s ft value (s) 


S'u (%) 


5ns (%) 






5c (%) 






nucleu 


Damgaard 


SM-SW 


SM-HF 


RHF-RPA 


RH-RPA 


IVMR 


T* = 

io c 


-1 : 

3041.7(43) 


1.679(4) 


-0.345(35) 


0.046 


0.175 


0.225 


0.082 


0.150 


0.008 


14 


3042.3(11) 


1.543(8) 


-0.245(50) 


0.111 


0.330 


0.310 


0.114 


0.197 


0.015 


22 Mg 


3052.0(70) 


1.466(17) 


-0.225(20) 


0.153 


0.380 


0.260 






0.031 


34 Ar 


3052.7(82) 


1.412(35) 


-0.180(15) 


0.285 


0.665 


0.540 


0.268 


0.376 


0.064 


T z — i 

26 Al 


: 

3036.9(9) 


1.478(20) 


0.005(20) 


0.182 


0.310 


0.440 


0.139 


0.198 


0.041 


34 C1 


3049.4(11) 


1.443(32) 


-0.085(15) 


0.326 


0.650 


0.695 


0.234 


0.307 


0.064 


38 K 


3051.9(5) 


1.440(39) 


-0.100(15) 


0.370 


0.655 


0.745 


0.278 


0.371 


0.077 


42 Sc 


3047.6(12) 


1.453(47) 


0.035(20) 


0.414 


0.665 


0.640 


0.333 


0.448 


0.091 


46y 


3049.5(8) 


1.445(54) 


-0.035(10) 


0.524 


0.620 


0.600 






0.106 


bu Mn 


3048.4(7) 


1.444(62) 


-0.040(10) 


0.550 


0.655 


0.620 






0.122 


b4 Co 


3050.8(10) 


1.443(71) 


-0.035(10) 


0.613 


0.770 


0.685 


0.319 


0.393 


0.139 


B2 Ga 


3074.1(11) 


1.459(87) 


-0.045(20) 


1.339 


1.48 


1.21 






0.175 


V4 Rb 


3084.9(77) 


1.50(12) 


-0.075(30) 


1.422 


1.63 


1.42 


1.088 


1.258 


0.235 


X 2 /n d 


(statistical experimental uncertainties only) 
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as a function of the mass number A. As an example, his 
expression in the microscopic model was 



5n = i8.o x iq- 7 a 5/s 



(22) 



We calculated the values of 5c from this equation for 
the thirteen best known superallowed transitions. The 
results, which are listed in column ten of Table HI were 
used as part of our comparative tests of all models. 



V. TEST RESULTS 

We have now set the stage for applying the test de- 
scribed in Sect. IIIII As explained there, our procedure 
for each of the six models is to compare that model's set 
of calculated 5c values (listed in Table HJ) with the set of 
values obtained from Eq. ([5]) and, using the method of 
least squares with Ft as the adjustable parameter, to op- 
timize the agreement between them. In effect, the 5c val- 
ues from Eq. ([5]) can be thought of as the "experimental" 



values: they incorporate the experimental ft values from 
Ref. [lj, as well as the small calculated correction terms, 
S' R and 5ns from Ref. j3| (also listed in Table HJ . The 
parameter, Ft, is a normalizer that allows each model 
to be tested for its success in obtaining a constant Ft 
value (i.e. in agreement with CVC), without regard for 
whether or not that Ft value ultimately satisfies CKM 
unitarity. The normalized x 2 f° r each least-squares fit - 
expressed as x 2 l n d, where rid is the number of degrees of 
freedom - thus yields a figure of merit for the model used, 
with smaller x 2 l n d values indicating better agreement. 

Although we take the measured ft values from Ref. [l| , 
strictly speaking the /t-value uncertainties quoted in 
that reference do not correspond to normal distribu- 
tions. Each ft value has three experimental inputs - 
energy, half-life and branching ratio - and each of these 
inputs typically includes a number of measurements of 
that quantity. The survey authors adopted the proce- 
dures used by the Particle Data Group [13| and, for any 
cases in which the measurements when averaged yielded 
a normalized x 2 greater than one, they increased the un- 



certainty on the average by a scale factor equal to the 
square root of the normalized % 2 ■ This conservative ap- 
proach leads to uncertainties on the ft and Ft values that 
are larger than would be the case for purely statistical re- 
sults. For our present purposes we have set all the scale 
factors in Ref. [l| equal to 1 and obtained new uncer- 
tainties on the ft values, which are normally distributed 
(at least to the extent that the uncertainties assigned by 
the authors of the original measurements were predomi- 
nantly statistical) . It is these re-determined uncertainties 
that appear in the second column of Table HI 

Obviously the uncertainties assigned to the theoretical 
radiative corrections 5' R , 5ns aud $c (if any) are not nor- 
mally distributed statistical quantities. Therefore, in our 
first least-squares test, we used only the re-determined 
uncertainties for the ft values and no uncertainties at all 
for any of the theoretical terms. The results for x 2 / n d 
appear in the first row below the main body of Table HI la- 
beled "statistical experimental uncertainties only" . Since 
this analysis uses only normally distributed uncertain- 
ties, we can proceed to evaluate a confidence level for 
each model. 

We follow the Particle Data Group [13[ in defining the 
confidence level (or "p- value" ) as being 



CL 



Pn d (X 2 )d X 2 , 



(23) 



where P nd {x 2 ) is the % 2 probability distribution function 
for rid degrees of freedom, and Xo is the value of x 2 ob- 
tained for a particular hypothesis - in our case, for a 
particular isospin-symmetry-breaking model. With this 
definition, the confidence level represents the probability 
that the x 2 f° r a valid hypothesis could exceed the value, 
Xo; actually obtained for the specific hypothesis being 
tested. More loosely, in our application the confidence 
level quoted for a particular model can be interpreted as 
the probability of that model being a valid one; i.e. of 
it being consistent with CVC. We express each CL as a 
percent on the next line in the table. 

We then present the results of a second least-squares 
analysis, in which we included uncertainties on the the- 
oretical radiative corrections, 5' R and 6ns', of course we 
retained the rc-dctcrmincd uncertainties already incor- 
porated for the ft values. The resulting x 2 l n d values 
appear in the second-to-last line in the table. These re- 
sults are also illustrated in Fig. [1] 

Finally, we list the results from a third least-squares 
analysis. Although two of the models, SM-SW and SM- 
HF, include theoretical uncertainties on 5c values in their 
original publications, the other four do not; so, to test 
them all on an equal footing, we have not used uncer- 
tainties on any of the model calculations in our first two 
analyses. We consider this to be the fairest approach. 
However, we have also examined what happens to our 
intercomparison if all the calculated 5c values are as- 
signed the same uncertainties as those originally quoted 
for the SM-SW calculations |3|. The values for x 2 / n d 



resulting from this third analysis appear in the bottom 
line of Table [D 

The most obvious outcome of these analyses is that 
only one model, SM-SW, produces satisfactory agree- 
ment with CVC, having x 2 / n d = 1.2 and CL = 26% 
in the properly statistical analysis. All of the other 5 
models have confidence levels well below 0.5%. Because 
the two other analyses included non-statistical uncertain- 
ties on the theoretical correction terms in addition to the 
statistical experimental ones, their values of x 2 l n d are 
substantially lower, but the relative ranking of the six 
models is approximately preserved: in all cases the SM- 
SW model is by far the best. It is remarkable that the 
model which becomes second best when the theoretical 
uncertainties are included is the earliest and arguably the 
most primitive one. Its success evidently stems from its 
treatment of the radial mismatch between the parent and 
daughter states, which accounts rather well for the sharp 
increase in 5c between Z = 12 and Z = 16, and between 
Z = 26 and Z = 30. It is perhaps equally striking that 
the most recent IVMR model fails to reproduce the trend 
of the data or any of its characteristic features. 



VI. CONCLUSIONS 

Evidently the shell-model with Saxon- Woods radial 
wave functions, SM-SW, is the only model tested 
that yields isospin-symmetry-breaking corrections which, 
when combined with the experimental ft values, produce 
Tt values that agree with the CVC hypothesis over the 
full range of Z values. This, of course, does not prove 
that the SM-SW model is correct in every way; however, 
it docs demonstrate that the other models in their present 
form cannot be used to extract a number for V uc i and to 
test CKM unitarity. As we note in Sect. [H] if the Tt val- 
ues are not consistent with one another, then their aver- 
age has no defined significance since either the symmetry- 
breaking model is wrong or CVC itself has failed. 

There is a second model, SM-HF, which has many 
promising features. As can be appreciated from an ex- 
amination of Fig. [TJ its relatively large x 2 is due to its 
failure to match the experimental 5c values for the cases 
with Z > 30. If we were to restrict ourselves only to 
the lighter cases, then the model would agree well with 
CVC. This difference at the highest Z values between the 
SM-SW and SM-HF model calculations has been known 
for 15 years, having first been pointed out by Ormand 
and Brown [|| even before the decays of the highest- i? 
emitters, 62 Ga and 74 Rb, had yet been precisely mea- 
sured. Prompted by the results reported here, we are 
currently examining whether this feature of the SM-HF 
model (as described in Sec. IIVC[) is sensitive to the par- 
ticular Skyrme interaction used [13]. We have by now 
sampled 12 different interactions and have also added a 
pairing term to the interaction, turning the calculation 
into a Hartree-Fock-Bogolyubov one. However, under no 
circumstances have we been able to produce agreement 



1.75 


i i i i 
Damgaard 


1 1 


^-v 1.25 
5 0.75 




A 


0.25 


j jIm- 


■Tfi 



1.75 
^ 1.25 
oo 0.75 

0.25 



SM-HF 


4 


► 



1.75 


RH-RPA (DD-ME2) 

« 


► 


^-v 1.25 




- 


to° 0.75 




- 


0.25 




- 




Y\ , * , 1 , , , , 



5 10 15 20 25 30 35 

Z of daughter 



1.75 


SM-SW 


J " 


^ 1.25 




/ 


to 0.75 


i 


► /§ • 


•-t-fi 


0.25 




- 



1.75 



1.25 



to" 0.75 
0.25 



RHF-RPA 


(PKOl) t - 








5 10 15 20 25 30 35 

Z of daughter 



FIG. 1: Isospin-symmetry-breaking correction, 8c, in percent units plotted as a function of atomic number, Z, of the daughter 
nucleus. The solid circular points with error bars are the values of 8c obtained from Eq. (JSJ), with the experimental ft values 
and the values of S' R and Sns (and their uncertainties) all taken from Table [I] In effect, we treat these as the "experimental" 
Sc values. The X's joined by lines represent the 8c values calculated by the various models described in the text and identified 
in the upper left of each graph. The value of Tt in Eq. ([5]) has been adjusted in each case by least-squares fitting to optimize 
the agreement between the "experimental" 8c values and the calculated ones. The corresponding values of \ 2 l n d are listed in 
the second-to-last row of Table |U 



with experiment over the full range of Z values. It is 
important to realize that both the SM-SW and SM-HF 
models use identical spectroscopic input, so it would ap- 
pear that the high-Z discrepancy is inherent to the SM- 
HF model itself. However, it must be admitted that too 
little spectroscopic information is known in this region to 
fully characterize the required model space. Calculations 
with larger model spaces and improved Hamiltonians are 
certainly to be encouraged. 

Fortunately, it is the successful SM-SW model that 
has principally been employed to calculate the 5c values 
used in the most recent data survey [l[. As was argued 
in that survey, the consistency of the Tt values was a 
powerful validation of those calculated correction terms 
and justified the subsequent derivation of V u d- However, 



in actually deriving V u d and its uncertainty we incorpo- 
rated the SM-HF calculations as well, even though we 
knew that model had a much poorer x 2 - Our rationale 
was one of conservatism. We enlarged the uncertainty 
assigned to the average J-t value to cover both sets of 6c 
calculations in order to be safe by including some pro- 
vision for systematic theoretical uncertainties. Whether 
we continue this practice in future is not yet decided. 

For now, though, we know that there are, as yet, no 
comparably successful competitive models. More impor- 
tant, we also have a protocol for testing future models, 
which is evidently very sensitive to the validity of the 
model. Furthermore, even though it is only the relative 
Z-depcndcnt variations in 8c that are being tested, it 
would surely require a pathological fault indeed in the 



theory to allow the observed nucleus-to-nucleus varia- 
tions in 5c to be reproduced in such detail while failing 
to obtain the absolute values to comparable precision. 

With this perspective, it is now informative to con- 
sider the points raised recently by Miller and Schwenk 
[10J . who claim that the SM-SW model is based on a for- 
mally incorrect interpretation of the isospin ladder oper- 
ator. They claim that this "incorrect" usage must have 
led to incomplete results for 6c, but they do not produce 
any "exact" calculations with which to compare. Instead, 
they identify a term involving radial excitations, which 
they consider to be missing from the SM-SW model, and 
proceed to evaluate this term under simplifying assump- 
tions. They assume that the radial excitations are domi- 
nated by mixing with states having one more radial node 
(c/. the Damgaard model) and, further, that the relevant 
excitations arc dominated by the isovector monopolc res- 
onance (cf. the IVMR model). Under these conditions 
they find that this "missing" term almost completely can- 
cels the SM-SW-model result and, although they produce 
no numbers, they state that this would result in 6c val- 
ues comparable in magnitude to the IVMR-model results, 
or even smaller. Clearly such a result would disagree at 
least as strongly with CVC as does the IVMR model. 
Therefore if any term is really missing from the SM-SW- 
model calculations, the test results presented here show 
that it must either be independent of Z or else very small; 
otherwise the data would become inconsistent with the 
CVC hypothesis. Considering that the Coulomb force is 



the principal source of isospin symmetry breaking, it is 
highly unlikely that any large component of 5c could be 
Z independent. 

From an experimental point of view, the results in 
Fig. Q] clearly demonstrate the importance of precisely 
measured ft- values. For example, the very precise val- 
ues for 26 Al m (plotted at Z = 12, the atomic number 
of its daughter) and 34 C1 (see Z = 16) contribute very 
significantly to the overall \ 2 f° r each model fit. Equally 
important though are the ft values for transitions that 
exhibit large values for 5c- The most obvious examples 
are the decays of 62 Ga (sec Z = 30) and 74 Zr (see Z 
= 36): their 5c values differ enormously from those for 
the transitions with Z<54 and this difference plays an 
important role in differentiating one symmetry-breaking 
model from another. More measurements of both types 
would be much welcomed in this context. 
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